Magnetic noise around metallic microstructures 
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We compute the local spectrum of the magnetic field near a metallic microstructure at finite temperature. Our 
main focus is on deviations from a plane-layered geometry for which we review the main properties. Arbitrary 
geometries are handled with the help of numerical calculations based on surface integral equations. The mag- 
netic noise shows a significant polarization anisotropy above flat wires with finite lateral width, in stark contrast 
to an infinitely wide wire. Within the limits of a two-dimensional setting, our results provide accurate estimates 
for loss and dephasing rates in so-called 'atom chip traps' based on metallic wires. A simple approximation 
based on the incoherent summation of local current elements gives qualitative agreement with the numerics, but 
fails to describe current correlations among neighboring objects. 

PACS numbers: 44.40.+a Thermal radiation -05. 40. -a Fluctuation phenomena, random processes, noise, and Brownian motion 



I. INTRODUCTION 

Thermal motion of charge carriers in a metallic object cre- 
ates a randomly fluctuating magnetic field in the object's 
vicinity. These fields are relevant for many applications like 
high-precision measurements of biomagnetic signals', nu- 
clear magnetic resonance microscopy^, and miniaturized traps 
for ultra cold atoms^ "*. The original purpose of Purcell's influ- 
ential 1946 paper^ was to point out that these magnetic fields 
have a spectral density that by far exceeds the Planck law for 
blackbody radiation at low frequencies. In fact, if only black- 
body fields were present, magnetic dipole transitions between 
atomic or nuclear levels would never happen on laboratory 
time scales. The near fields sustained by material objects (that 
play the roles of sources and cavity) give the dominant con- 
tribution. Indeed, these near fields contain non-propagating 
(evanescent) components that are thermally excited as well 
and that dominate over free space radiation''-^. Phrased in an- 
other way, the dipole transition rate is enhanced because the 
elementary excitations in the metal provide additional decay 
channels2i2. 

We focus in this paper on accurate calculations of mag- 
netic field noise that are able to describe objects of arbi- 
trary shape. Such objects occur, for example, in magnetic 
microtraps where complex networks of metallic wires create 
electromagnetic potentials with typical scales in the micron 
ranged. The behaviour of the field spectrum, as the metallic 
geometry is changed, is far from intuitive. The spectral density 
increases with the material volume for small structures, but 
this trend saturates as soon as the typical scale gets larger than 
the penetration length (skin depth) of the fields in the material. 
It has even been found that a thin metallic layer can produce 
less noise than a half-space, depending on the ratios between 
observation distance, layer thickness, and skin depth^'^Mi. 
Experiments in the field of biomagnetism have shown signifi- 
cant changes when a metallic film is cut, at constant volume, 
into stripes'. Calculations for these cases necessarily require 
numerical methods to describe the propagation of magnetic 
fields both in vacuum and inside a metallic structure. This is 
the main topic of this paper. We also discuss previously de- 
veloped approximations for planar structures and within the 
magnetostatic regime where analytical calculations are possi- 



ble. Our numerical methods are restricted here to two spatial 
dimensions (2D) where efficient solutions of Maxwell equa- 
tions can be found with the help of boundary integral equa- 
tion&i^iiiiiii^. Our approach can also be combined with any 
other numerical method for field computations, permitting to 
cover three-dimensional cases as well. 



The results we find can be summarized as follows. A planar 
structure (infinite lateral size) creates equal magnetic noise for 
all components of the magnetic field vector. This does not ap- 
ply in three dimensions, but is specific to the two-dimensional 
setting we focus on here. Finite metallic objects show a strong 
anisotropy: the noise occurs preferentially along 'azimuthal' 
directions circling around the object. Increasing the amount 
of metallic material does not always give larger noise, in par- 
ticular when the geometrical size becomes comparable to the 
skin depth. We find qualitative agreement with measurements 
of Ref, 1 where thermal field fluctuations are reduced when a 
metallic object is spUt into disconnected pieces. The surface 
impedance approximation, that provides an accurate descrip- 
tion of metallic reflectors for far-field radiation, is shown to 
be not reliable for observation distances shorter than the skin 
depth. Finally, qualitative (albeit not quantitative) agreement 
is obtained between our numerical data and an approxima- 
tion based on the incoherent summation of fields generated by 
thermal current elements filling up the metallic volume. This 
method has been used in the interpretation of previous exper- 
iments^^'^^. Our results are, to our knowledge, the first quan- 
titative test of this approximation in a nontrivial geometry. 



The paper is organized as follows. We first review the link 
between the thermal radiation spectrum and classical dipole 
radiation (Sec|II]i. Planar structures are analyzed in Sec Hill 
using the angular spectrum representation. We demonstrate 
in particular the isotropy of the magnetic noise spectrum and 
discuss the accuracy of the surface impedance approximation. 
SeclTVlis devoted to our numerical scheme and to the results 
for single and multiple objects of rectangular shape. 
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II. MAGNETIC NEAR FIELD NOISE 

A. Local noise power 

The fluctuations of the thermal magnetic field B(r,t) are 
characterized by their local power spectral density (the Fourier 
transform of the autocorrelation function) 



By(r;c^) = /dTe"^-(B,(r,t)B,(r,i + r)). 



(1) 



Higher moments are not needed for our purposes, and the av- 
erage field (at frequency ui) vanishes as is typical for thermal 
radiation. We assume the field to be statistically stationary 
so that the spectrum ([T]l does not depend on t. The diago- 
nal tensor components Bii{r; uj) give the spectrum for a given 
cartesian component Bi{r) or polarization direction. Previ- 
ous work has shown a strong dependence on the position r 
near a metallic microstructure: power laws being the typi- 
cal behaviour in the frequency range where the wavelength 
A — 2ttc/uj is much larger than the typical distances. For the 
temperature dependence, see Eq.© below. Our parameters 
of interest are: normal metallic conductors with temperatures 
above a few K, r in the micron range and A of the order of 
centimeters or larger (w/27r < 10 GHz). This upper limit on 
frequency corresponds to the strong magnetic dipole transi- 
tions in typical alkali atoms.) In this regime, the frequency 
dependence of the noise spectrum is weak and occurs via the 
material response (permittivity e{uj)). A characteristic length 
scale is the field penetration length (skin depth) 6 defined by 



27r 

T 



(2) 



where a{u!) is the conductivity and /io the vacuum perme- 
ability. The second expression is based on the Hagen-Rubens 
approximation Ime(u;) w a{uj) / {eqCli) 3> |Ree((jj)|. Within 
the Drude model for a metal, this is verified at frequencies 
much smaller than the charge carrier relaxation rate (in the 
10^^ range at room temperature). For highly conducting 
materials (Au, Ag, Cu), this results in a skin depth of the order 
of 100 ^m(w/27rMHz)-i/2^ 

To compute the magnetic correlation spectrum, we use the 
fluctuation-dissipation theorem which is valid at thermal equi- 
librium (temperature T)'^ 



2h 



Im^„(r,r;tj) 



(3) 



where Qi.j{Y,Y' \lo) is the Green function for the magnetic 
field, i.e., the field generated at r by a point magnetic dipole 
located at r' and oscillating at the frequency lo, Bi{r,t) ~ 
Qij (r, r'; u!)^j e~'"* + c.c. This is actually a familiar result: 
the imaginary part of trtj(r, r; gives the local density of 
magnetic field modes, and the temperature dependent prefac- 
tor in Eq.(l3]l their average occupation number The basic ben- 
efit of this formula is that it holds also for the full correlation 
tensor and even near material objects that absorb the field or 
generate thermal radiation. Generalizations to the nonequilib- 
rium case exist (fields produced by a 'hot object' surrounded 



by a 'cold' environment)'^ '^, but are not needed for our pur- 
poses (see remarks in Sec|V]l. We also note that the temper- 
ature dependent prefactor in Eq.© reduces to 2kBT/uj for 
T > 0.1 K. In this limit, the order of field operators in the 
correlation function ([T]l becomes irrelevant. (The order we 
have adopted yields the rate of a magnetic dipole transition 
i — > / with energy difference Ef — Ei — huj.) 



B. Magnetic dipole radiation 

We are thus led to solve the following electrodynamic prob- 
lem: find the complex magnetic field amplitude B(r;aj|/x) 
created by a monochromatic point dipole fi{t) — fie' 
c.c. located at position r'. We then compute 



Gy(r,r';w) 



dBi{r;uj\fi) 



(4) 



In the limit r ^ r' this field becomes the singular 'self field' 
and requires a cutoff in wavevector space. Its imaginary part 
is cutoff-independent, however, and given by ImB(r';Lj) = 
^QLu^fj,/ {Gttc^) (in three-dimensional free space). 

The field B = B(r;Ci;|/x) can be found from the vector 
potential A that solves the inhomogeneous Maxwell equation 



V X V X A - fcoe(r)A = ^qV x fiS{r - r'), 



(5) 



where fco = ^/c The right-hand side is the current density 
corresponding to the magnetic dipole. There is no free charge 
density and we work in the gauge E = icjA. 

We now focus on the following geometry (figlU right): the 
position r' of the source (i.e., where the magnetic noise spec- 
trum is actually needed) is located in vacuum, and the metal- 
lic microstructures are filling a domain T) where Im e(r; uj) — 
a{v; ut) /{squj) is nonzero (and large). 



B(r') 

^ ©: vacuum 
( B : metal 



|X\ B(r') 



\ \ j(r) 



v : vacuum 



'D : metal 
a(r) 



FIG. 1: Sketch of the considered geometry: (left) current fluctua- 
tions in a microstructure generate magnetic field fluctuations B(r') 
at a position position r' outside it. (right) The magnetic noise spec- 
trum is calculated from the magnetic field radiated by a point mag- 
netic dipole ^ located at r'. T> and V': domains where the conduc- 
tivity a(r; lu) is nonzero or zero, respectively. 

The outside domain is called V . There, the vector po- 
tential satisfies an inhomogeneous Helmholtz equation with 
wavenumber fco. All length scales we consider (distance 
dipole-microstructure d, object size) are much shorter than 
the wavelength so that fcg is actually very small and can be 
neglected in a first approximation. This is the magnetostatic 
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regime. (The finite value of ko is, of course, at the origin of 
the nonzero magnetic LDOS in free space.) We cannot make 
the magnetostatic approximation in V because there, we have 
a wavenumber fco-\/e(r;^) — (1 + ™d the (local) 

skin depth S{r) is one of the characteristic length scales at 
hand. The fields in the domains T) and V are connected by 
the usual matching conditions: the components of A tangen- 
tial to the boundary are continuous, and B is continuous (the 
material is non-magnetic). 

Eq.(|5]l provides a unique solution subject to the boundary 
condition that at infinity, the field behaves like an outgoing 
wave. In three [two] dimensions, this corresponds to a vec- 
tor potential proportional to e'^^^^'/s [e^''"'' / y/s] in the free 
space domain V' when the distance s = |r — r'| oo to 
the source becomes large compared to A. In the magnetostatic 
limit fco — * 0, the free space asymptotics is actually never 
reached at finite distances. The relevant boundary condition is 
then the same as for the scalar potential of an electric dipole: 
the vector potential goes to zero like [like 1/s] in three 
[two] dimensions, respectively. 

Since we deal with a metallic object with |e| ^ 1, it 
is tempting to perform the calculation based on the surface 
impedance boundary condition. The latter links the tangen- 
tial components of magnetic field and vector potential by 
Bt = -iLoZA^, where loZ = (1 + i)/S. Note that this is 
a local relation that can only hold if the scale of variation of 
the fields on the object surface is much larger than the skin 
depth 6. In the present study, a point-like source illuminates 
the object with its near field [Abuik(r ^ r') in Eq.(|7l)], and 
this field shows a typical extension of the order of the object- 
source distance d. The surface impedance approximation is 
hence expected to break down for d ^ 6. We shall confirm 
this explicitly for the planar structures discussed in the follow- 
ing SecHni 



III. RESULTS: LAYER 

A. Two dimensions 

In this paper, we focus on a two-dimensional (2D) geome- 
try to simplify the numerical calculations described in SeclIVI 
The magnetic moment is chosen in the computational plane 
(the xy-plane), as shown in FiglU Adapting the wave equa- 
tion (|5]l to two dimensions, we find that the vector potential 
has a single nonzero component that points out of the plane. 
We then work with a scalar function A{r) = A{x, y) that 
solves 

V^A + kle{r)A = (My^^' - t^xdy') 5{r - r') (6) 

In a homogeneous medium ('bulk'), the solution with the ap- 
propriate boundary conditions is 



^ ifJ^yda;' - fixdy') Ho{ko^/e\r~r'\) (7) 



Auik(r-r') 

where Hq is the Bessel function of the third kind (Hankel 
function), usually denoted -ffo^"* = Jo + i^o- i^Froni this. 



we get the magnetic field by taking the 'curl', = dyA, 
By — —dxA. The resulting self field in free space is 

ImB(rV) - iMofc^M (8) 

provided the dipole /i. is real. In the magnetostatic limit, this 
field is negligibly small. The bulk solution (|7]l then goes over 
into 

A (r- r-M- A^o {x - x')ny ~ {y ~ y')^^ 



2tt 



r — r' 



This equation describes the field with which the dipole 'il- 
luminates' the sample. Note that it is scale-free: the typical 
'spot size' on the microstructure is only determined by the 
distance d between dipole and top surface. 



B. Reflected fleld 

In this section, we consider that the boundary of the 
medium is the plane y = 0; the field at the source point 
r' = (0, d) is then related to the (Fresnel) reflection coeffi- 
cients from the surface. We expand the solution to Hq.® in 
plane waves (wavevector k parallel to the boundary) and have 
above the medium (y > 0): 



+ 00 



A{x,y) = tltlilJ-xdy' ~ IJLydx') 



dk e''^^^-^') 
27r 2k 



(^e-''\y-y'\ + r{k)e-''^y+y"^^ (lo) 



where k = ^Jk"^ — k^. (The square root is chosen such that 
Re K > 0.) The coefficient r{k) describes the reflection of the 
field from the medium boundary. It is given by the Fresnel 
formula 



T(k) — ?'half space (fc) 



- ^k^ - 2i/<52 

(11) 

for a medium with skin depth 6 filling the half-space y < 0. 
For a layer (thickness h) on top of a substrate, we have 



^laycr 

(fc) 



^top ~t~ T'bottom^ 



-2K„,h 



1 



^top ^bottom 6 



(12) 



where rtop = ^haif space is given by Eq.O and rbottom de- 
scribes the reflection from the layer-substrate interface. It is 
given by Eq.lfTTT) with the replacements n ^ Km, Km ^ Ks = 
{k^ — Es^o)^^^ where Ss is the substrate permittivity. 

All the relevant information for the magnetic noise power 
is contained in the reflection coefficient r{k). In fact, when 
the integral in ( fTOl ) is performed and the imaginary part taken, 
it turns out that the reflected waves (second term) dominate 
over the free space contribution (first term) by at least a fac- 
tor A^J/d^ ^ 1. This is connected to the fact that the rel- 
evant wavenumbers k for our problem are of the order of 
l/{y + y') = l/(2c?) which is much larger than fco- We can 
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hence apply the approximation k « |fc|. The reflection coef- 
ficient ( fTTT l for the metalHc half-space then depends only on 
the parameter kS. For the metallic layer geometry, we fo- 
cus for simplicity on a substrate whose conductivity is much 
smaller than in the metal. The influence of the substrate has 
been studied in Reffiol: already a ratio of 10 to 100 between 
the substrate and layer conductivities is sufficient to make the 
substrate behave like vacuum. We then have rbottom ~ —''top 
in Eq.jnil. 



C. Polarization dependence 




0.001 



100 



Let us analyze first the dependence on the orientation of 
the source dipole. If //, is perpendicular to the medium (only 
Hy 0), the reflected field is given by 



dfc P 
2n2K 



r{k)e 



ik(x-x') -K{y+y') 



(13) 



FIG. 2: Reflection coefficients dill 112b for thin and thick metallic 
layers. We plot the imaginary part only. The dashed lines represent 
the formulas of Table|l] The wavenumber is scaled to the inverse skin 
depth 1/5. For the thick layer, h — oo. We take the conductivity of 
gold at room temperature and a frequency ui/2tt ~ 1.1 MHz in all 
figures. This leads to the value S — 71fim and a vacuum wavelength 
\/5^3.S X 10^ 



The Umit r — ^ r' yields an imaginary part 

— — r(fc)e-2-'^. (14) 

Repeating the calculation for a parallel dipole, we find for 
lmBx{v'\fix) the same expression as Eq.(fT4l). and conse- 
quently the noise spectrum is isotropic, Bxx — Byy This is a 
remarkable property of a laterally infinite structure in 2D. (In 
3D, the polarization perpendicular to a planar interface has a 
noise power twice as large as the parallel polarizationJi^.) We 
show below that a significant polarization anisotropy arises 
above a metallic wire of finite width. 





k < 1/S 


l/S < k 




Im rhalf space (fc) 


kS 


1 

2^252 






k < 1/S 




1/h < k 


Im naycr(fc) 


-l + {l + i)kS 
™ l + ifc(5Vft 


h 
kS^ 


1 

2fc2(52 



TABLE I: Asymptotic approximations to the reflection coefficients 
from a half-space and a layer. We distinguish between thin (thickness 
h <^ 5) and thick layers (h > S, 'half space'). The first and second 
columns (thin layer) overlap in an intermediate fc-range (see Fig|2}. 
The magnetostatic limit fco >C fc is taken throughout. These formulas 
are plotted as dashed lines in Fig[2] 



D. Wavevector dependence 

The reflection coefficient ( fTSl l is plotted in Fig|2]for typi- 
cal layers. Consider first a thickness larger than a few skin 
depths. One observes a maximum value of its imaginary part 
(relevant for the magnetic LDOS) when the decay constant 
K K k is matched to l/S. This is confirmed by an asymptotic 
analysis whose results are given in Table I] (See Refs|20||2ll 
for details on the asymptotic expansion.) One of the two lim- 
iting cases (namely fcg <C fc <C 1/S) corresponds precisely 
to the surface impedance approximation where the reflection 
coefficient (fTTT ) is approximated by 

r(fc")«-l + (l + i)|fc|5 (15) 

Here, the skin depth is much smaller than the lateral period 
and the field barely penetrates into the material. Fig|2] and 
Table U show strong deviations in the opposite regime k ^ 
1/6 that is relevant at distances d ^ S. 

Consider now a layer much thinner than the skin depth. 
^^From Fig|2l different regimes can be read off that are sep- 
arated on the fc-axis by the scales h/S'^ <^ 1/h, as can be 
seen in Fig|2] It is worth noting that for small fc, thin layers 



show even larger losses [Imr(fc)] than thick ones; the maxi- 
mum is shifted towards the smaller value fc ^ h/S^ and has 
a larger amplitude. This behaviour has been recognized be- 
fore in magnetic noise studies in the kHz range^. In the in- 
frared range, it is also well known that the absorption by a 
metallic layer can be optimized at a specific thickness. (See, 
e.g., Ref 22 for incident far-field radiation where |fc| < fcg.) 
Conversely, for a given thickness h and dipole distance d, the 
magnetic noise power shows a maximum as the skin depth is 
changed^ '". This 'worst case' occurs when the characteristic 
wavevector l/d is matched \.o h/5'^. 

E. Distance dependence 

The asymptotics in fc-space translate into power laws for 
the dependence of the magnetic power spectrum Bii{d;uj) 
on distance d, as shown in Figl3] In fact, the integrand in 
Eq.(fT4li peaks around k ~ l/{2d), and the result of the in- 
tegration is determined, to leading order, by the behaviour 
of r{k) in this range. We thus find the power laws summa- 
rized in Table inland visible in Fig[3] We use as convenient 



5 



03 

I 1 
Q. 

CD 

0.01 

o 

o 10-4 

E 

o 
c 



10" 



\ 


1 


thick 


1 


thin 







0.01 0.1 1 10 100 
distance [skin depth] 

FIG. 3: Local magnetic noise power Bii{d;uj) vs. distance from 
metallic layer in double logarithmic scale (2D calculation). Top 
curve [blue]: thick layer; bottom curve [red]: thin layer. The dashed 
lines give the leading order power laws of Table Ull The thick curves 
arise from the numerical integration of Eq.dlOb. the thin curves are 
an interpolation formula described in the text. 
The magnetic noise power is isotropic above a planar structure in 
2D (the perpendicular and parallel field components have the same 
power). It is scaled to ^oksT / {ujS'^), and the distance is scaled to 
the skin depth 5. Thin [thick] layer: h = 0.01 5 [3 5]. 
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S^/h < d 
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log(/i/d) 


h 




2-K 


2^ 


4TThd^ 



TABLE II: Power laws for the magnetic noise spectrum in two di- 
mensions above a half space and a thin metallic layer (dashed lines 
in Fig|3j. The noise spectrum is given in units of ixqUbT / {ujS^). 



gions between the power laws. 

Finally, the surface impedance approximation gives a mag- 
netic noise that is represented in Figl3] by the dashed line 
close to the 'thick layer' for d > 5. The agreement with 
the full calculation in this range is expected: the 'illuminat- 
ing field' is getting more and more uniform on the scale of the 
skin depth. At shorter distances, the surface impedance ap- 
proximation severely overestimates the noise level because it 
cannot describe properly the field variations on scales smaller 
than 5. For the thin layer, the conventional surface impedance 
approach gives a wrong result even \f d > 5 because top 
and bottom surfaces do not decouple from each other. This 
can be repaired using effective (thickness-dependent) surface 
impedances, see, e.g.. Refill and citations therein. 



IV. RESULTS: FINITE SIZE OBJECTS 

We now describe numerical calculations that we have per- 
formed to assess the importance of the finite lateral size of the 
metallic structure. This is particularly relevant, for example, 
in atom chips where a continuous metallic layer is etched to 
define wires that can be addressed with different cuiTents^ii24. 
It is actually desirable to minimize the amount of metallic ma- 
terial, leaving just a few wires to create the fields for atom 
trapping. In fact, it has been argued that the magnetic noise 
power roughly scales with the metallic volume as long as the 
characteristic distances are smaller than the skin depth^i^l. For 
laterally finite structures, this claim as well as other calcu- 
lations have been based so far on approximate methods that 
fail to reproduce even the planar layer to within a factor of 
two or threeii*2i2^i22. The numerical results we describe here 
are a first step towards an accurate estimate of magnetic noise 
power near structures of finite size. 



unit in all the plots the noise level fjiQk-BT / {ojS'^). Normal- 
ized to blackbody radiation (in 2D free space), this level is 
(2/(fco(5))2 1.5 X 10^^ at LI MHz for gold at room temper- 
ature, a striking illustration of the Purcell effect^. A common 
trend is that the magnetic noise power increases as the metal- 
lic medium is approached. As the distance d is getting much 
smaller than the thickness h, thin and thick layers behave the 
same, as expected. At larger distances, but still smaller than 
the skin depth, the noise power is proportional to the volume 
of metallic material, hence to the layer thickness^i^^. This 



trend is reversed for d > Sy^d/2h where thin layers give a 
larger noise level than thick ones^iiS. 

A reasonably accurate approximation that interpolates be- 
tween these power laws can be found by performing the k- 
integral using the asymptotic formulas of Table U in their re- 
spective domains of validity. The result is a sum of incomplete 
gamma functions r(n, x, x') (n = 0, 1, 2) whose arguments 
are, for example, x ~ dh/6'^, x' ~ d/5, or d/h (thin lines in 
Fig|3]and AppendixlAll. We have checked that the asymptotics 
of the gamma function reproduce the power laws summarized 
in Table There are regimes where the sub-leading terms 
give significant corrections, in particular in the transition re- 



A. Boundary integral equations 

Within the assumption of near field radiation being in equi- 
librium with the metallic object, we compute the noise power 
from the magnetic Green function in Eq.©. The magnetic 
field radiated by a point source and reflected by the object 
solves the wave equation (|6]l. We reformulate the wave equa- 
tion in terms of boundary integral equations. This has been 
described elsewhere' ^''^''^-'^ , and we quote only the basic for- 
mulas here. Our unknowns are the nonzero component A of 
the vector potential and its normal derivative F = dA/dn = 
n • VA on the object surface S, where n is the outward unit 
normal to S. Both quantities are continuous (actually, F is 
equal to the tangential magnetic field) and can be found from 
the system of integral equations 



A{t) = 2Abuik(r - r') 



2 ^da(x) 
A(r) = 
2 ^da(x) 



dGx 



Gi(r-x)F(x)-^(r-x)A(x) 



Gar-x)F(x)-^(r-xM(x) 



(16) 



(17) 
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where Abuik(r — r') is given in Eq.©. Both the observation 
point r and the integration points x are taken on the object 
boundary S here, da(x)) being the surface element. We use 
the scalar Green functions [see Eq.©] 



Ge(r) 



-Ho{koV^\r\). 



(18) 



If we would take the magnetostatic limit, Gi(r) 
— (27r)^^ log |r|, the Green functions in vacuum and in the 
medium would differ (in sub-leading order) by a constant, 
leading to inconsistencies. We avoid this by retaining the fi- 
nite value of fco even for the vacuum Green function. The 
integrals ( fT6l [TtT i are to be understood as principal values to 
handle the singularities of the Green functions as x ^ r. We 
discretize them on a finite element decomposition of the object 
boundary, as described in Ref 12. The resulting linear system 
is solved with standard numerical tools. Once the fields A, 
F are known on the surface, the reflected field at the source 
position (r') can be found from 



— ^da{x) 



Gi(r'-x)^^(x) 



dGi 
dn 



(r' - x)A(x) 



(19) 



Note that Gi and dCi / dn are both essentially real here (the 
imaginary parts scale with ko). The magnetic noise, via 
Im B(r'), is thus determined by the imaginary parts of A and 
F on the object boundary. This is not surprising since the 
induced current density is aE = iujaA. 



B. Single wire 

We have solved the integral equations for rectangular wires 
of thickness h and width w. In a first step, we have validated 
our numerical scheme by comparing flat, wide wires {w 3> h) 
to the infinite layer results of Sec|lII] Typical plots are shown 
in FigH] where the magnetic noise power (symbols) is plot- 
ted vs. the distance d above the wire centre. Good agree- 
ment with the analytical results for an infinitely wide wire 
(solid lines) is only obtained at short distance, where for ge- 
ometrical reasons the wire appears wider. At distances above 
20 ixm, the deviations start to grow. In all the plots, we take 
a skin depth (5 = 70 /im. The slow convergence in the limit 
w ^ oo can be attributed to the long-range behaviour of the 
fields; this is more pronounced in two dimensions compared 
to three. Note in particular the strong splitting between the 
two polarization directions for the thick wire that does not oc- 
cur above an infinitely wide wire in 2D (Sec lIIICl l. Interest- 
ingly, the ?/-component shows more noise above a thick wire 
while this tendency is reversed above a thin wire. This polar- 
ization anisotropy could provide a tool to improve the lifetime 
in a magnetic trap: one orients the static trapping field paral- 
lel to direction of the strongest noise. (In fact, trap loss and 
spin flips are induced by magnetic fields perpendicular to the 
static trap field.) The choice of a trapping field along the weak 
noise direction is favorable if one wants to reduce the dephas- 
ing rate of the trapped spin states (generated by fluctuations 
of the Larmor frequency, see Ref.^). 




10 -20 30 40 50 60 70 

distance [nm] 

FIG. 4: Magnetic noise spectra Ba (d) vs. distance d above the centre 
of thin and thick metallic wires. Symbols + (x): numerical calcula- 
tion for the Byy (Bxx) component (see Fig[T). Solid lines: infinitely 
wide wire (layer), as computed in Sec|Ill] Dashed lines: incoherent 
summation (thin layer only, upper curve: Bxx), see Sec. lIVCl Thin 
wire: width and thickness 200 x 7 fim; thick wire: 200 x 160 /im. 
The skin depth is 5 = 70 fim. 



Another finite-size effect is shown in Fig|5]where the posi- 
tion is varied parallel to the top surface of a thin wire. Above 
the centre of wide wires, the noise levels are nearly constant 
(not shown). Beyond the wire edges, one observes a sharp 
drop in Bxx, with a characteristic scale fixed by the distance. 
The //-component shows a broad maximum near the edge that 
is more pronounced for naiTow wires. This is due to a gradu- 
ally changing direction of maximum noise that is 'azimuthal' 
with respect to the object, as expected for magnetic fields gen- 
erated by currents flowing perpendicular to the computational 
plane (see inset of FiglUl. We find the direction of maximum 
noise by looking for the eigenvectors of the 2x2 matrix Bij . 
This matrix can be shown to be symmetric (reciprocity), and 
the eigenvectors that are not orthogonal very near to the wire's 
corner, are an artefact of our numerical method that converges 
very slowly at these points. 

In Figs|6] |2l the thickness of the wire is changed with the 
observation point remaining above the centre. We observe 
an approximately linear increase with the width that saturates 
slowly. We also note that Bxx (left) levels off faster than Byy 
(right). The difference between Fig|6]and Fig|7]is the distance 
of observation: at short distance (Fig|6]l, the largest widths 
show a noise power fairly close to the planar layer limit (cf. 
the symbols at the right end). At distances comparable to the 
skin depth (Fig|7]), the deviations from the planar layer limit 
(symbols) are still large. Note also that the noise has dropped 
in amplitude and that the increase with width is slower. 



C. Incoherent sununation 

This behaviour can be qualitatively understood using the 
'incoherent summation' approximation developed in Refl23l 
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FIG. 5: Magnetic noise spectra Bii{d) vs. lateral position, at a fixed 
distance d. The arrows mark the edges of the wires. Symbols x [+]: 
spectrum Bxx [Byy] parallel [perpendicular] to the top face of the 
wire. Skin depth: 5 = 70 nm. 

Left panel: thickness and width are 20 x 200 /im (wide wire) and 
20 x 20 (narrow wire). Distance d = 10 /im. Right panel: thick- 
ness and width are 1 x 20 /xm, distance d — 3 fim (see dotted line of 
inset). Inset: illustration of anisotropic noise near the wire edge. The 
crosses are oriented along the polarization vectors that show max- 
imum and minimum noise, the 'arm lengths' being proportional to 
the rms noise. The magnetic field noise is dominantly azimuthal, 
with field lines circling around the wire. The dotted line (d — 3 /im) 
shows the positions scanned through in the right panel. 




FIG. 6: Magnetic noise spectra vs. the width of a rectangular wire, 
(left) x-polarization, parallel to the top face; (right) j/-polarization. 
Symbols: numerical calculations; solid lines: incoherent summation 
approximation (see Sec lIVCl ). The symbols on the right margin give 
the values for an infinitely wide wire (layer). The observation point 
is located above the wire centre, at a distance d = 10 nm. The 
wire thickness is 7 /im (thin) and 160 nm (thick). Skin depth: S — 
70 /im. 
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FIG. 7: Same as Fig|6] but at an observation distance d = 75 fim. 
Results from the incoherent summation are not shown, as they 
strongly deviate. 



FIG. 8: Noise power generated by three wires, as a function of dis- 
tance (see inserted sketch, with the dashed line illustrating the ob- 
servation points). The wires have a quadratic cross section 20 /im x 
20 fim and are separated by a gap of 20 /im. The noise is measured 
above the center of the central wire. Left panel: horizontal polar- 
ization, right panel: vertical polarization. Symbols: numerical re- 
sult; solid line: incoherent summation. For comparison is shown: 
a single wire of same cross section (dashed line) and a wide wire 
20 fim X 80 fim with approximately the same volume (dash-dotted 
line). Skin depth 5 — 70 fim. 



the metallic volume is broken into mutually incoherent point 
current elements whose magnetic fields are computed within 
magnetostatics and neglecting the presence of the metallic ob- 
ject. We give the resulting formulas for two dimensions in 
Appendix|Bj the integrals are solved by special functions for 
a wire with rectangular cross section. The solid lines in Fig|6] 
demonstrate that incoherent summation gives a reliable ap- 
proximation if the skin depth is the largest length scale (not 
true for the thick wire). The noise power always increases 
with the metallic volume within this approximation, however, 
and it may also happen that a wider wire produces a slightly 
weaker noise (x-polarized curve for a thin wire in Fig|6ll. This 
is qualitatively similar to the trend of Fig[3]where a thick layer 
can produce less noise than a thin one at distances larger than 
the skin depth. The polarization anisotropy is also qualita- 
tively reproduced by the incoherent summation method, al- 
though Bxx is overestimated. In fact, due to damping on the 
scale of S, not the entire volume of the thick layer contributes 
to the noise. The dashed lines in Fig|4]and further calculations 
show that the quantitative agreement is systematically better 
for the field component perpendicular to the nearest metal sur- 
face (here, Byy). 



D. Multiple wires 



This is the generic situation in miniaturized magnetic traps 
('atom chips') with wires being defined by etchings in a metal- 
lic layer We consider three wires of identical cross section 
and smaller than the skin depth. We show in Fig[8]the depen- 
dence on the vertical distance, above the central wire. Our 
results interpolate smoothly between a single narrow wire 
(d ^ w) and a single wide wire (d ^ 2w), as could have 
been expected. In fact, the three geometries give nearly the 
same noise in the azimuthal (Bx) polarization. The incoher- 
ent summation approximation overestimates this noise com- 
ponent (similar to Fig|6]). We attribute this to correlations be- 
tween the current fluctuations in neighboring wires that are not 
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FIG. 9: Same as FigjS] but at fixed distance d = 10 fim, scanning 
thie lateral position, a; = is above the center of the central wire). 



captured by incoherent summation. On the other hand, this 
approximation gives an excellent agreement for the weaker 
noise component By. 

When we shift the observation point laterally, along the axis 
connecting the wire centres, we get Fig|9l The stronger B^- 
polarization shows maxima of noise above each wire, as ex- 
pected. In the Sy-polarization, a maximum occurs in the gap 
between the wires. This conforms to the general trend of 'az- 
imuthal noise' illustrated in Fig |5] (inset). It is also interesting 
that above the central wire (x — 0), three wires generate less 
noise than only one and also less than a single wide wire (ap- 
proximately a merger of the three). This observation goes into 
the same direction as the experiments reported by Nenonen 
and co-workers' where a reduction of thermal magnetic fields 
was achieved by cutting a metallic film into stripes. We at- 
tribute this behaviour to negative correlations between the cur- 
rents in neighboring wires brought about by the propagation 
of the magnetic field between them. In fact, the noise could 
only increase if the wire currents were strictly uncorrected. 

The performance of the incoherent summation approxima- 
tion (solid lines) can be clearly seen, the trends being similar 
to FiglD good agreement for the _Bj^ -polarization, overestima- 
tion of the perpendicular case due to the neglect of correlation 
effects. 



for the discrepancy between experiment and theory observed 
in Reffiol although our method (restricted to 2D) do not per- 
mit quantitative predictions of trap lifetimes. 

We now comment on possible extensions of this analysis. 
Our framework is also able to provide an approximate de- 
scription of superconducting structures. In fact, since we deal 
with magnetic field fluctuations at a finite frequency, there is 
always some penetration into the superconductor or, equiv- 
alently, a finite resistivity. This can be attributed to a frac- 
tion of carriers in a normally conducting state. Calculations 
for superconductors of planar geometry have been reported 
in Refs il0ll29l with applications for miniaturized atom traps 
in mind. More accurate descriptions require one to solve the 
London equations at finite frequency inside the superconduc- 
tor, using for example the two-fluid model^. 

We recall that we use in this paper a local version of Ohm's 
law. For very pure metallic films, the ballistic transport of 
charge carriers implies a nonlocal respons o-^''^^'^^'^^ . This may 
be particularly relevant for wires defined by doping a semi- 
conductor, but would require major changes for the numerical 
approach of Sec II VI 

Finally, a brief remark on non-equilibrium settings. Con- 
sider an isolated metallic object held at a temperature different 
from its surroundings (materialized by the vacuum chamber 
walls, for example). By applying the generalized Kirchhoff 
relations (see, e.g., Refil7b. the radiation arriving at the ob- 
servation point can be split in two parts: one is proportional 
to the product of the power a test dipole emits into the far 
field and the temperature of the surroundings; the second part 
is proportional to the dipole radiation power absorbed by the 
object and the metal temperature. At the sub-wavelength dis- 
tances of interest for this paper, one can show that the second 
part is dominant and that the error made in using the same 
temperature for metal and surroundings is small. The equilib- 
rium calculation we have focused on here is then sufficient. 



V. CONCLUSIONS 

We have described in this paper numerical and analytical 
results for the thermal fields surrounding a two-dimensional 
metallic object of arbitrary cross-section. The role of the skin 
depth (5 as a characteristic length scale has been highlighted. 
At distances smaller than S, the spectral noise power roughly 
scales with the volume of the metallic material (Figs.|3]|6]|3- 
We have reviewed a simple method based on this idea, the 'in- 
coherent summation approximation' . It systematically overes- 
timates the noise power in one of the two field polarizations, 
but otherwise reproduces the main features as long as the skin 
depth is the largest scale. The strong polarization anisotropy 
that we have found suggests strategies to minimize loss or de- 
coherence due to thermal magnetic fields, as observed in re- 
cent experiment&i^iZ^i^: this can be done by suitably choos- 
ing the direction of the static trapping fields. We have also 
shown that the noise power can be significantly non-additive 
when dealing with multiple objects. This could be relevant 
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APPENDIX A: UNIFORM APPROXIMATION 

The integral over k in Eq.(fT4ll can be performed analyti- 
cally if the power-law approximations of Table J] are used for 
the reflection coefficients. We split the integration range at the 
crossing points between the power laws and sum the contribu- 
tions. The full expression is cumbersume, and we quote here 
only the most complicated case, the thin layer in the range 
fco <C fc ^ 1/(5 (first column of Table [III. The integral can be 
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handled with the formula 

fe2 



dfc 



(-1 + {l + i)kS)kc 
1 + ikt 



-2kd 



in the xy-plnne that we compute in the magnetostatic approx- 
imation and ignoring the presence of the embedding metal. 
The latter point is the key approximation made. This gives a 
magnetic noise spectrum (integrated over a unit length along 
z) of the order of 



(l-i)^r(2,2fcid,2fc2d) 



(Al) 



2d 



T{l,2kid,2k2d) 



+ y e-2'''/*r(0, 2kid - 2id/t, 2k2d - 2id/t) 

where t — S'^/h and / = (1 + i)^'^/^ + i/^- I" this case, 
ki — ko and fc2 = {tS)~^^'^. The incomplete gamma function 
is defined by 



T{n,b,a) = T{n,b) — T{n,a) 

oo 

T{n,a) = [dtt"-^-* 



(A2) 
(A3) 



Using the asymptotics of the gamma function-'^, we get the 
power laws in Table HIl The logarithmic behaviour arises from 



a<l: r(0,a) w -loga. 



APPENDIX B: INCOHERENT SUMMATION 



(A4) 



We outline here the adaptation of the incoherent summation 
method of Ref 23 to two dimensions. The thermal spectrum of 
the current density is given at low frequencies (huj ^ ksT) 
by 

(j*(x;cj)j(x';w')) = 27r5(tj - tj')2fcBrCT(x; cj)5(x - x') 

(Bl) 

This spectrum is already integrated over a unit length in the 
z-direction (parallel to the current) along which the current 
density is assumed to be uniform (two-dimensional geome- 
try). In this formulation, a is (the real part of) the 3D con- 
ductivity that we assume local, as reflected by the spatial S- 
correlation. We only take into account currents parallel to the 
z-direction. Each current element generates a magnetic field 



Se 



with cross correlations given by 



Sr 
TT 



dx^dx^ ■. TTj — — 



(B2) 

(B3) 
(B4) 



where V is the volume occupied by the metal. The 'geomet- 
rical tensor' Yij is dimensionless (a specific 2D property) and 
depends only on the ratio of observation distance and object 
size. It does not involve the skin depth, of course. 

For a microstructure with rectangular cross section, an ob- 
server located above the center of the structure sees a noise 
power 



Sb 



1 

2^ 



arctan(- 



1 

'2^ 



Byy{d;uj) 



1 

2^ 



IX 

ImLi2( -) 

y-y . 



arctanf^ — 
x' 



1 

'2^ 



ImLi2( t) 

y-y . 



y' = -h 


(B5) 



y' = -h 




J y' = -h 



where Li„(-) is the polylogarithm and we have used the nota- 
tion 
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fia,c)-fia,d)-fib,c) + fib,d) 

(B6) 
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